Dust Control in Finite Air Volumes at Zero Gravity - Mean-Field Like Analysis 
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We study a simple ID model of dust rods, with mean size \i, passing through a parallel ID 
alignment of pores as a problem of clogging of a filter by dust grains; fi is kept less than the pore 
size, s. We assume that the filter is "sticky", characterized by some parameter < A < 1, which 
means that dust grains slightly smaller in size than s can get trapped in the pores. Our analyses 
suggest that the number of clogged pores, iV c j, grows in time as t^° s cx where v = v((j,, A) is a 
non-universal exponent that depends upon the dust size distribution and filter properties. 

PACS numbers: 47.55.Kf,47.55.Mh,89.20.Bb 



Common household dust typically comprises of partic- 
ulate matter such as carbon particles, skin flakes, dust 
mites, fiber, hair and other irregular shaped objects with 
length-scales between 10 _2 //m and 10 2 ^m (see, for ex- 
ample, ,1]). Typically, the larger sized hairs/fibers act 
as nucleation centers for agglomeration of fine scale dust 
and the eventual development of the familiar dust bun- 
nies. Presence of dust bunnies reduces the effectiveness of 
dust filters. We study the problem of control of dust at all 
length-scales to ensure that dust bunny growth becomes 
difficult. In particular, we consider dust control at zero 
gravity in an air volume that is finite and cannot be re- 
plenished. Our mean-field like analysis is inspired by the 
problem of air quality maintenance in the International 
Space Station. The present study has the potential to 
facilitate the design of filters which can retain air-quality 
below specified thresholds, for long periods. 

Significant work has been done to describe the pres- 
sure drop in fluids with fine particulates that is flowing 
through filters with given specifications (see, for exam- 
ple, |3 )■ The analyses have been performed typically via 
rate equations using state variables, akin in spirit to ther- 
modynamical analysis. The resultant blocking equations 
reveal the behavior of pressure drop with time as filtrate 
volumes are modified and filter criteria are tuned [2j. 
Limited information appears to be available on specifics 
such as time-dependent degradation of filters as functions 
of dust size distributions and filter pore distributions. 

The present work focuses on filter clogging as a func- 
tion of time. We refrain from addressing the dust ag- 
glomeration process in this Letter. Fine scale dust is 
continuously generated in most inhabited environments 
and naturally agglomerates. Dust mites (~ 10 2 /mi), 
plentiful in most environments, tend to feed on smaller 
length scale dust. Presence of dust mites therefore leads 
to a natural integration of smaller length scale dust into 
larger length scale dust. 

In our analyses, we will be considering dust size distri- 
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butions such that the mean dust size is comparable to the 
filter pore size. This turns out to be the regime of interest 
in which one can understand how filter and dust proper- 
ties lead to the time evolution of the clogging process of 
the filter. Our studies show that dust size distributions 
can affect the details of the clogging process. 

The Model System: The present study is a II?, mean- 
field-like, analysis. Dust grains are irregular in shape 
and come in all sizes. In our ID model, we only consider 
the size of the dust grains, i.e. effectively consider rod 
shapes. Consideration of all sizes in ID is accomplished 
by assuming a Gaussian distribution of dust sizes. Our 
studies suggest that an exponential tail to the distribu- 
tion or a Lorentzian distribution will affect the rate at 
which filters clog. However, the major qualitative fea- 
tures remain the same for the different distributions con- 
sidered and, therefore, Gaussian distribution will form 
the basis for our study. 

Randomly sized rods from the distribution, assumed 
to be aligned vertically for easier geometric conceptual- 
ization, are allowed to impinge at random positions on 
a vertical ID filter with vertical slits or "pores". The 
walls between the pores are considered to be of negligi- 
ble width. Rods of all sizes travel at fixed speeds towards 
the filter. We assume that our time steps are sufficiently 
fine and/or the air is sufficiently dilute in dust content 
such that a single rod reaching the filter is a reasonable 
approximation. This justifies measurement of time, t, in 
terms of particle release events at 'unit' time intervals of 
At; n th particle is released at t = n x At. Not all parti- 
cles contribute to clogging. We will only count clogging 
events, and measure clogging rates, with respect to com- 
plete clogging of a pore. The number of clogged pores 
at any time is denoted as N c i. One would expect N c i to 
initially grow rapidly in time. But as time progresses, it 
will become progressively difficult for the rods to find the 
next available pore to clog. Thus, the temporal growth of 
N c i will eventually saturate. To best describe the man- 
ner in which N c i grows, we define a new parameter to 
characterize the instant at which the N th clogging event 
happens and call it t c ^ og . This also leads us to the param- 
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N — t N _^), which is a waiting time between 
1)"* and N th clogging. As time progresses and 
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FIG. 1: In Fig. ^a), A is kept fixed while varying fi and, in Fig. UJ[b), vice versa. It is clear from these graphs that t c N ° 9 
oc Nd"^'*' and this is explicitly shown in Fig. Q^c), where all the graphs have been collapsed onto a single curve through 
v = fj,- 62 x (1 - A)- 2 . 



the filter gets progressively clogged, we expect the aver- 
age value of t to grow in magnitude. 

The mean of rod size, r, is referred to as fi and is 
kept smaller than the filter pore size, s. The standard 
deviation of rod sizes is taken to be unity. It is to be un- 
derstood, however, that the variance will affect filter life, 
filtering efficiency etc. and will be discussed later. We 
study cases of filters with fixed pore sizes, as well as filters 
with a Gaussian ditsribution of pore sizes; filter pore size 
is initially kept fixed in the results reported below, and is 
to be considered so unless otherwise mentioned. Ideally, 
when r < s, rods go through open pores, and, when r 
> s, they block open pores while, also, partially blocking 
adjacent pores by the amount by which r > s. However, 
we allow a rod with r slightly < s to stick to a pore; this 
models the various sticking tendencies induced by diverse 
physical mechanisms, and also generates the well docu- 
mented 'bridge formation' 0], where several particles link 
together to arch over a pore. For this, we define a stick- 
ing threshold, < A < 1, which prescribes the rod length 
that can block an open pore as a fraction of s. If A = 1, 
only rods with r > s can block a pore, and hence this is 
the 'non-stick' limit. If A < 1, particles with r < As will 
pass through the pore, without sticking. Typically, 0.5 < 
A < 0.9 in our calculations, [i is set arbitrarily equal to 5 
units, and s to 9 units. 

We also allow for the formation of "cakes", another 
commonly observed phenomenon (see, for example, @|), 
with the additional deposition of rods on top of previ- 
ously deposited rods, which may not contribute to ad- 
ditional blocking unless they extend over adjacent open 
pores. Rods are not assumed to hit the barrier between 
pores with any significant probability and get trapped by 
the same; this eliminates the need to consider the partial 
blocking of adjacent pores by a rod that drapes over the 
barrier. Further, adjacent pores are considered only to 
one side (in the direction of increasing array index) for 
convenience; these effects do not compromise the valid- 
ity of the final results. In the following, we present our 
results on the clogging process of the ID filter. It may 
be noted that a 2D filter can be thought of as an array 



of adjacent ID filters. In the limit of rods being incident 
at the same rate on each ID 'unit' filter, the results of 
our present study are expected to remain applicable in 
the 2D case. 

Results: We consider a filter with Np pores at time 
t = 0; Np, in the following, is kept at 10 5 for adequate 
statistical robustness. In Fig. ^ we plot the time taken 
for the N th clogging event, t c ^° 9 , versus the number of 
clogged pores, N c i, for various choices of the two control 
parameters, \i and A, in log- log mode. In Fig. ^a), A 
is held fixed while varying /x and, in Fig. M ^ s held 
fixed while A is varied. All the cases are characterized by 
two different regimes of power law growth where, firstly, 
there is a smooth increase in t c ^ 9 with N c i, followed by an 
abrupt changeover to a region of infinite slope; the latter 
marks the divergent times, on the average, required to 
clog the final few pores, which are scattered randomly on 
the ID filter. Clearly, the interesting physics is in the 
first region, and we see that t c ^° 3 oc N c i v , where v = v(fJ,, 
A). The scaling pattern is made explicit in Fig. [Jc), 
where we have collapsed all the various graphs for the 
different /x's and A's onto a single graph, and it emerges 
that 

A) =M Q1 x (1.0- A) Q2 . (1) 

where a.\ = —0.62 and ct2 — —0.2 for this case of Gaus- 
sian distribution of rod sizes impinging on a filter with 
fixed pore size. We found that these exponents are some- 
what "fragile" and can vary depending on the nature of 
the rod and filter pore size distributions. For example, 
with an exponential distribution of rod sizes and a fil- 
ter of fixed pore size, the exponents are a± — —0.2 and 
a 2 = —0.1. Nevertheless, the scaling behavior is robust. 

Figs.^a)-(&) reveal that the time required to clog all 
the pores tend to grow dramatically when N c i becomes 
sufficiently large. We recall that the time period between 
successive cloggings is r. In Fig- 121 we present data of the 
distribution of t, P(t), versus r, for two typical cases; the 
plots are in log-log mode and only r's upto 10 4 are shown. 
In one case, P(t) decays smoothly while, in the other 
case, the decay is preceded by a plateau region where the 
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10° x 10 4 

FIG. 2: A plot of the distribution of waiting times, P(t), 
against r, in a log-log plot, shows a plateau region preceding 
the decay region when fi is reduced; the graphs are with fixed 
pore size, s = 9.0. A similar plateau region results if A is 
reduced, keeping a fixed. This indicates the enhanced role for 
smaller particles in the distribution, which are characterized 
by smaller r's, through the incremental clogging mechanism. 



r's have the same density of occurrence. In general, it was 
observed that the plateau region increases with relative 
decrease of /x with respect to s. Also, the same effect is 
generated if A is decreased. In both these cases, we note 
that the smaller particles have an enhanced role to play 
in the clogging process and the plateau region is increased 
due to the phenomenon that we call incremental clogging 
(IC). 

One can see that the role of the big (r > As) rods is nec- 
essary for the initiation of the clogging process; note that 
these come from the tail of the distribution, have smaller 
probabilities and, hence, larger r's associated with them. 
Once the process has been initiated, the smaller rods can 
contribute to additional clogging through the IC process, 
i.e. in the next stage, only a rod of size bigger than the A 
fraction of the remaining open portion of the pore is re- 
quired to further block the pore and so on. These smaller 
rods have larger probabilities and smaller r's associated 
with them. Also, if, at any stage, the rod that sticks on 
is big enough to protrude into the adjacent pore, it can 
reduce the size of that pore, which initiates IC in this 
adjacent pore. While the big rods contribute to clogging 
when IC is active, their role is especially significant in 
the initial stages of the clogging process. IC increases 
the proportion of small r's and, as a corollary, when the 
role of smaller particles is enhanced either by reduction in 
/x or by reduction in A, IC becomes a more active mecha- 
nism. It is to be emphasized here that the random nature 
of the clogging events mandates non-zero probability of 
obtaining small r's even in the final stages of the clogging 
process; indeed, this has been clearly seen in plots of r's 
versus N c i Q. 

At the macroscopic level, it is of interest to probe the 



average rate of clogging at the time of clogging of the 
N th pore. We define this average rate of clogging via K c i 
= Nd/t^ 09 . This quantity is shown in Fig. |3| plotted 
against N c i, for various choices of control parameters, /x 
and A. In Fig. [21 a), A is held fixed and [i is varied while, 
in Fig. |2I b), /x is held fixed and A is varied. The data 
suggests two different power-law growth regimes of K c i 
as a function of N c i . There is a smooth increase in K c i in 
the early stages followed by a steep decrease in the rate 
when there are only a few pores left to be clogged; the 
former regime corresponds to the functional period in the 
life of the filter. It turns out that the plots in Fig. ^a) 
and Fig. can be unified, as shown in Fig. c) . The 
scaling relation used on the data in Fig.|3Ja) and Fig.|!Jfr) 
turns out to be 

K cl = k x N*- v . (2) 

where k = exp(/(/x, A)), and /(/x, A) = —16.3 + 1.04/x + 
4.66A. The relationship with respect to N c i is to be ex- 
pected given the definition of K c i and Eq. (1). However, 
further complexity of the scaling relationship via the pref- 
actor k results because we are considering a rate quantity 
in a problem with discontinuous kinetics. 

It is easy to calculate (see, for example, Q), based 
on the probabilities for particles of size r > ko~, where 
a is the standard deviation of rod size distribution and 
k any positive value, the probability of finding a rod of 
size r > s as against finding a rod of size r > As; it 
turns out that the latter is greater by three orders of 
magnitude if A = 0.8, with /x = 5.0, a = 1.0 and fixed 
pore size s = 9.0 j4(. Differences in probability translate 
to comparable differences in t c ^ 9 . Clearly, variance of the 
distribution also controls the rate of clogging. Similar 
differences in probability can be achieved by varying /x 
with respect to s. 

If we introduce a distribution of pore sizes around a 
mean value, instead of having a fixed pore size, t c ^ 9 de- 
creases; in this comparison, the same value of the fixed 
pore size was assigned for the mean. If we introduce a 
second mode at a lower value in a corresponding uni- 
modal Gaussian distribution, t c ^° 9 increases because of 
the enhanced proportion of smaller sized rods; in this 
comparison, a value of 3.0 was used for the second mode 
with the first mode at 5.0. An exponential distribution, 
which has a higher proportion of bigger rods than a Gaus- 
sian, leads to faster clogging of the pores. These results 
are summarized in Fig. 0] in a log-log plot of t c ^° 9 ver- 
sus N c i. A corresponding ordering, with respect to the 
distributions, can be observed in plots of K c i versus N c i 
plots; note that larger times translate to lower K c i values 
in this case. 

Concluding Remarks: Our studies indicate that ad- 
justing the filter pore size (or mean of the same in case 
of a pore size distribution) to be slightly higher than the 
mean of the dust particle distribution is an important 
step in obtaining efficient filtering as well as longer filter 
lifetimes. The variance as well as the type of distribu- 
tion of the filter and dust particles are also important 
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FIG. 3: Plots of K cl (= N cl /t%° 9 ; cf. Fig. □ versus N c i, in log-log mode; in Fig. ^j[a), A is kept fixed while varying fi and, in 
Fig.|^j6), vice versa. A scaling relationship akin to Fig.^c) may be suspected; this is indeed so, as shown by Fig.^Jc), where all 
the graphs have been collapsed onto a single curve through proper scaling, with k = exp(/(/i, A)), f(u, A) = — 16.3+1.04/X+4.66A. 




FIG. 4: Effect, on clogging rates, of different rod size dis- 
tributions, with and without filter pore size distribution are 
compared. The notation in the legend is: BM- Bimodal Gaus- 
sian (two modes being at 5.0 and 3.0), G- Gaussian (mean 
at 5.0), E- Exponential (mean at 5.0); FP- fixed pore size (s 
= 9.0), PD- distribution of pore sizes with mean at 9.0. It 
is seen that faster clogging takes place with a pore size dis- 
tribution in general, and bimodal distribution is the slowest 
and exponential distribution fastest in clogging the filter. 



determinants. The scaling law that we propose in the 
relationship of t c ^ 9 and N c i, with respect to the con- 
trol parameters fi and A, can play an important role in 
suitably choosing the filter parameters. Even though the 
scaling exponent v is somewhat "fragile" with respect to 
different distributions of dust particles, the fact that its 
magnitude is only altered by a factor of the order of 2, 
rather than by an order of magnitude, in going from a 
Gaussian distribution to an exponential distribution, is 
very encouraging. A key parameter for further analy- 
sis is the resistance to the flow caused by filter clogging 
which requires modelling of fluid flow; the relationship of 
cakes to the model parameters would also be interesting 
to study. The present study is expected to be appli- 
cable to any problem where small particles can lead to 
clogging of any fluid flow. The model may have broad 
applicability with respect to problems where the concept 
of incremental clogging can be applied. 
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